Wave-packet dynamics at the mobility edge in two- and three-dimensional systems 
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We study the time evolution of wave packets at the mobility edge of disordered non-interacting 
electrons in two and three spatial dimensions. The results of numerical calculations are found to 
agree with the predictions of scaling theory. In particular, we find that the return probability 
P{r = 0,t) scales like t~^^^'', with the generalized dimension of the participation ratio D2. For 
long times and short distances the probability density of the wave packet shows power law scaling 
P{r,t) oc f-02/d^D2-d^ rpj^g numerical calculations were performed on network models defined by 
a unitary time evolution operator providing an efficient model for the study of the wave packet 
dynamics. 
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The time evolution of wave packets is intimately re- 
lated to the transport properties of an electron system. 
If a wave packet does not spread unboundedly for long 
times the system is an insulator. On the other hand, if a 
wave packet spreads diffusively for long times the system 
has a finite conductivity characteristic of a metal. In a 
seminal paper Anderson observed that disorder can lead 
to the absence of diffusioiad|-| With the advent of the scal- 
ing theory of localizationnm it was realized that the ex- 
istence of this Anderson transition between localized and 
metallic states depends essentially on the dimensionality 
and the symmetries of the system only. In one dimen- 
sion any amount of disorder leads to localization of all 
states while in two dimensions for non-interacting elec- 
trons either r& strong magnetic field (as ia the quantum 
Hall effectETLl) or spin-orbit interactionErtHl are necessary 
for the existence of extended states. Three-dimensional 
system generally show an Anderson transition as function 
of disorder strength and Fermi energy. 

At the Anderson transition a system shows a behavior 
intermediate between localized and extended. A crucial 
assumption of the scaling theoryB is the finite value of 
the conductance G of the system at the transition. The 
conductance of a d-dimensional hyper-cube of linear di- 
mension L is related to the conductivity a by Ohm's law 
G = aL'^^^. The conductivity at the transition thus 
scales like L^^'^ and so does the diffusion coefficient D 
according to the Einstein relation a — e^pD, with the 
density of states p. The Anderson transition is not re- 
flected in thcpdensity of states p which is finite through 
the transitionQ. Accordingly, the diffusion coefficient is 
finite at a two-dimensional transition while there is no 
diffusive behavior in three dimensions. 

Apart from this scale dependence of the diffusion co- 
efficient a distinguishing feature of the Anderson tran- 



sition is the emergence of multi&j-actal correlations of 
the wavefunction amplitudeg3'Elrll3. The corresponding 
density-density correlator can be expressed in terms of 
a wave vector- and frequency-dependent diffusion coeffi- 
cient D{q,Lo), which in the DC limit becomes the diffu- 
sion coefficient discussed in the previous paragraph. For 
short distances and low frequencies, however, this diffu- 
sion coefficient refiects the sjbppg amplitude fluctuations 
within single eigenfunctionsBliJ. It should be noted that 
multifractal correlations are also present away from the 
Anderson txaasition on scales smaller than the localiza- 
tion lengthEjO. 

The scaling form of the diffusion coefficient-/i(c, w) has 
been discussed previously in the literaturelialllrEy. In this 
paper we study directly the time evolution of wave pack- 
ets. The physical information contained in both quanti- 
ties is the same. The advantage of our approach is that 
it obviates the need for a diagonalization of the Hamilto- 
nian. Instead the time evolution operator is applied re- 
peatedly. For systems defined by a Hamiltonian H this is 
not necessarily easier to do as the time evolution operator 
U = exp{—iHt/h) is an exponential of the Hamiltonian 
and the kinetic and potential terms in the Hamiltonian do 
not commute in general complicating the calculation of 
the exponential. To circumvent this problem decomposi- 
tions of the exponential haV|e-|been devisednil and applied 
to the problem at handcaiHj. On the other hand, for 
systems defined by a unitary operator that can be inter- 
preted as a time evolution operator the dynamics of wave 
packets is the most natural and convenient quantity to 
study. A class of systems that falls into this, category are 
the network moclck introduced by ShapirocJ and Chalker 
and Coddingtoncj. A network model is characterized 
by a unitary operator that describes the scattering of 
wavefunction amplitudes at the nodes of a network. Net- 
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work models have been constructed and studied for two- 
and three-dimensional systems with-jptiiogonal, unitary, 
symplectic, or chiral symmetrycZrEJM'E^I. The unitary 
operator of a network model can be interpreted as a dis- 
crete time evolution operator where a time steDjXprre- 
sponds to the typical time of a scattering eventEjS. 

The spreading of wave packets has also been discussed 
in connection with the nature of the spectrum of the sys- 
tem. Ketzmerick et al. found that the spreading of wave 
packets is determined by the multifractal dimensions D2 
and Z?2 of rthe local spectrum and the eigenfunctions, 
respectivelyEj. While in general these dimensions are in- 
dependent quantities characterizing a system, they are 
simply related in disordered system, D2 = dD2, since 
the global density of states is scale invariant and— con- 
nects length and energy scales via the dimension aEB. In 
this paper we will restrict our discussion to the case of 
the Anderson transition. 

After summarizing the scaling behavior of the prob- 
ability density of a wave packet, we present the results 
of numerical simulations for two- and three-dimensional 
network models and show that they are in accordance 
with the predictions of scaling theory. In particular, we 
calculate the return probability to the origin of the wave 
packet, and the shape of the wave packet for long times 
and short distances. We conclude the paper with a dis- 
cussion of the results. The scaling forms for the shape 
and moments of the wave packet as well as numerical 
results for the moments and jjxeliminary results for the 
shape have been presented inE3. 

We now summarize the predictions of scaling theory for 
characteristic features of a-wave packet. For the deriva- 
tion we refer the reader tcxB and references therein. 

Due to the occurence of multifractal eigenfunction fluc- 
tuations at the localization-delocalization transition, the 
shape of a wave packet differs strongly from the diffusive 
gaussian form. The long time, short distance behavior, 
r'' ^ i/?ifl-.with p the density of states, of P{r,t) is a 
power lawO, 

P{r,t) cct-^'/'^r"'-'^. (1) 

The exponent D2 is the generalized dimension of the in- 
verse participation ratio. For multifractal eigenstates it 
is smaller than the space dimension d. The term short 
distance has to be understood in the sense of r'^ <C t/fip. 
The distances on which Eq. (|^) holds are large compared 
to any microscopic scale and correspond for t — > 00 to 
the limit r ^ 00. For long times Eq. (|l|) thus describes 
the bulk of the wave packet up to an exponentially small 
tail. 

From eq. (|l|) it follows that the return probability to 
the starting point of the wave packet P{r = 0, t) decays 
with time as 

P(r = 0,i) oci-^^/'^. (2) 

This is the first quantity that we turn our attention 
to. It is the simplest quantity that reflects the multi- 
fractality of the eigenstates at the mobility edge. From 
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FIG. 1. Return probability for a wave packet at the quan- 
tum Hall critical point. The solid line is a best fit to the data 
with slope D2/2 — 0.76 ± 0.03. The dashed line shows the 
classical behavior p{t) oc t^^ . 

Fig. (|l|) we get the exponent D2/2 of the return prob- 
ability to be 0.76 ± 0.03, somewhat smaller but still in 
agreement with the value 0.81 ± 0.02 obtained by Huck- 
estein and Schweitzer pfor the quantum Hall transition in 
a tight-binding modeO. In the present calculation, 400 
realizations of the disorder were averaged for a system of 
100 X 100 scatterers. 

For a three-dimensional system the exponent D2/S of 
the critical power law differs even more strongly from 
the diffusive value of 3/2. The system that we stud- 
ied is a stack otxuiantum Hall networks coupled in the 
third dimensioiic3^EB. For vanishing inter-layer coupling 
the system shows the quantum Hall transition of the iso- 
lated two-dimensional systems, while for finite coupling 
the system shows a three-dimensional Anderson transi- 
tion with localized and non-localized regimes separated 
by mobility edges. 

The return probability for the three-dimensional net- 
work model at the Anderson transition is shown in 
Fig. (j^). For the time interval from 400 to 4900 time steps 
we obtain a power law with an exponent D2/3 = 0.43 ± 
0.04. For shorter times a crossover to two-dimensional 
behavior is observable in our data. Its origin is the 
strongly anisotropic character of our system of coupled 
two-dimensional networks. The coupling in the third di- 
mension is much weaker than in the planes so that three- 
dimensional spreading is observed only after sufhciently 
long times. The data presented in Fig. (||) were calcu- 
lated from 200 disorder realizations of a system of 50'^ 
scatterers. 

When comparing our result D2 — 1.3 ± 0.1 with val- 
ues published in the literature it should be kept in mind 
that our network model corresponds to a system with 
a magnetic field. For time reversal symmetric systems 
in the absence of magnetic fields somewhajLJarger values 
have been published previously: 1.7 ± 0.3t3, 1.45-1.8E3, 
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FIG. 2. Return probability for a wave packet at the An- 
derson transition in d = 3. The solid line is a best fit to the 
data with slope D2/3 = 0.43 ± 0.04. The dashed line shows 
the conventional behavior p{t) oc t^^''^. 

1.7±0.2@, 1.5±0.2ii. The values obtained previously in 
the absence of time reversal sy|iHmetry |do not differ sig.- 
nificantly from these: 1.7 ± 0.2^, « l.SEZi, 1.55 ± 0.150. 
Again our value is somewhat but not significantly smaller 
than these values. 

We now show that the bulk of a wave packet has a 
power law shape. For the two-dimensional quantum Hall 
system at the critical point the average probability den- 
sity of the wave packet is plotted in Fig. (^) for different 
times t. The different curves have been rescaled by a 
factor of t^^/^ in order to allow for a better comparison 
of densities at different times. The straight line shows 
the expected power law with exponent D2 — 2 = —0.48 
where the value D2 was taken from the fit to the return 
probability. The crossover from power law to exponential 
shape of the density distribution is seen to shift with in- 
creasing time to larger radii. 6 realizations of the disorder 
were averaged for a system of 300 x 300 scatterers. 

Figure (0) shows the shape of the wave packet in a 
three-dimensional system of 50^ scatterers. At the criti- 
cal point the power law with an exponent D2 — 3 = —1.7 
is observed after 20000 time steps (solid boxes). In the 
metallic regime the density is practically constant over 
the whole system after 2400 time steps (open boxes). In 
the localized regime exponential localization is observed 
after 20000 time steps (triangles). 

To summarize, we have studied the time evolution of 
wave packets at the mobility edge in disordered two- and 
three-dimensional electron systems. In order to charac- 
terize the wave packets we looked at the return proba- 
bility and the shape of the density distribution. Due to 
the scale invariance of the critical systems these quanti- 
ties show power law scaling. The exponents differ from 
the usual diffusive situation for two reasons. On the one 
hand, at the Anderson transition the conductance is scale 
invariant and not the conductivity or the diffusion coeffi- 
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FIG. 3. Average density of wave packets at the quantum 
Hall critical point. The different curves correspond to differ- 
ent times (1500, 3000, 6000, 9000, 12000, 15000 time-steps). 
The different curves have been rescaled by an amount pro- 
portional to t° ''^. 
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FIG. 4. Average density of wave packets 
in three-dimensional network models. Data are plotted: in 
the metallic regime after 2400 time steps (□), at the Ander- 
son transition after 20000 time steps (■), and in the localized 
regime after 20000 time steps (A). 

cient. As a result the spatial dimension d enters where in 
the diffusive case the exponent 2 of the diffusion equation 
enters. On the other hand, the short distance long time 
behavior is governed by multifractal density correlations 
within single eigenfunctions leading to the replacement 
of the spatial dimension d by the multifractal dimension 
Z?2 of the inverse participation ratio. 

We have confirmed the scaling laws and extracted the 
scaling exponents by numerical calculations. The calcu- 
lation were performed for network models of disordered 
systems that are defined by a unitary network operator 
that serves as a discrete time evolution operator. This 
provides a very convenient and efficient way to study the 



3 



time evolution of wave packets. 

The time dependence of the return probability to the 
origin shows power law scaling with the exponent —D2/d 
involving the multifractal exponent D2. Previously, this 
quantity was studipiin two- and three-dimensional tight- 
binding modelaljcJ'Lfl. We obtain values for D2 of 
1.52 ± 0.06 and 1.3 ± 0.1 in two and three dimensions, 
respectively. These values are somewhat smaller, espe- 
cially in three dimensions, than previously published val- 
ues. At present, it is not clear whether this difference re- 
flects genuinely different critical behavior in the systems 
or different strength of corrections to scaling. It is known 
that the two-dimensional network model has very small 
corrections to-^paling compared to the random Landau 
matrix modelMa. 

For long times t and short distances, r'^ -C t/fip, the 
shape of a wave packet becomes t-'D2/d^D2-d ^ ^■\^_ 
serve this behavior in the two- and three-dimensional 
network models. At larger distances the crossover to ex- 
ponential tails in thu-density distribution is observed. 

Ketzmerick et a/.Ej obtained under quite general con- 
ditions the expression P{r,t) (x t~^^r^^~'^ for the shape 
of wave packets, where D2 is the generalized dimension 
of the local density of states. These general results agree 
with our results for the Anderson transition since in the 
present case the generalized dimensions for the spectrum 
and the wave functions are not independent but pro- 
portional to each other, D2 = dD2cBtB. The origin of 
this simplification is the scale invariance of the global 
density of states that allows to relate energy and length 
scales via a single relevant energy dependent length scale 
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